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Spontaneous reconnection at a separator current layer. 

I. Nature of the reconnection 

J. E. H. Stevenson,^ and C. E. Parnell^ 

Abstract. Magnetic separators, which lie on the boundary between four topologically- 
distinct flux domains, are prime locations in three-dimensional magnetic fields for recon¬ 
nection, especially in the magnetosphere between the planetary and interplanetary mag¬ 
netic field and also in the solar atmosphere. Little is known about the details of sepa¬ 
rator reconnection and so the aim of this paper, which is the first of two, is to study 
the properties of magnetic reconnection at a single separator. Three-dimensional, resis¬ 
tive magnetohydrodynamic numerical experiments are run to study separator reconnec¬ 
tion starting from a magnetohydrostatic equilibrium which contains a twisted current 
layer along a single separator linking a pair of opposite-polarity null points. The result¬ 
ing reconnection occurs in two phases. The first is short involving rapid-reconnection in 
which the current at the separator is reduced by a factor of around 2.3. Most (75%) of 
the magnetic energy is converted during this phase, via Ohmic dissipation, directly into 
internal energy, with just 0.1% going into kinetic energy. During this phase the recon¬ 
nection occurs along most of the separator away from its ends (the nulls), but in an asym¬ 
metric manner which changes both spatially and temporally over time. The second phase 
is much longer and involves slow impulsive-bursty reconnection. Again Ohmic heating 
dominates over viscous damping. Here, the reconnection occurs in small localised bursts 
at random anywhere along the separator. 


1. Introduction 

Many highly energetic space physics processes involve 
magnetic reconnection, such as solar and stellar flares, coro¬ 
nal mass ejections, interactions between planetary and inter¬ 
planetary magnetic fields and substorms in magnetospheres, 
etc. Two key papers that explained the basics of generalised 
three-dimensional (3D) reconnection [Schindler et al., 1988; 
Hesse and Schindler^ 1988] provide the corner stones to all 
3D reconnection studies. In particular, one of their main re¬ 
sults reveals that, unlike 2D reconnection, 3D reconnection 
can either occur with or without null points. 

In 3D, reconnection has been shown to occur at a mul¬ 
titude of sites, such as at topological features including 3D 
null points (locations where all three components of the mag¬ 
netic field equal zero) [e.g., Craig et a/., 1995; Pontin et a/., 
2004, 2005; Pontin and Galsgaard, 2007; Priest and Pontin^ 
2009; Masson et aL, 2009; Pontin et al.^ 2011] and mag¬ 
netic separators (special field lines that link pairs of 3D null 
points) [e.g.. Priest and Titov, 1996; Longcope and Cowley, 
1996; Longeope, 2001; Haynes et al., 2007; Parnell et al., 
2010a, b; Wilmot-Smith and Hornig, 2011] or at the geo¬ 
metrical features known as quasi-separatrix layers (regions 
of magnetic field, which at one end are closely anchored, 
but at the other end are anchored far apart) [e.g., Priest 
and Demoulin, 1995; Demoulin et al, 1996, 1997; Aulanier 
et al, 2005, 2006] and in twisted or braided flux tubes [e.g., 
Calsgaard and Nordlund, 1996, 1997a, b; De Moortel and 
Calsgaard, 2006; De Moortel and Calsgaard, 2006; Wilmot- 
Smith and De Moortel, 2007; Browning et al, 2008; Hood 
et al, 2009; Bareford et al, 2013]. 
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The reconnection that we focus on here is that associated 
with magnetic separators. Such reconnection has been in¬ 
voked as an explanation for a number of specific observed 
solar flares [e.g., Longeope et al, 2005; Barnes et al, 2005], 
for flux emergence events [e.g., Parnell et al, 2010b; Mae- 
Taggart and Haynes, 2014] and in the heating of the quiet- 
Sun [e.g.. Close et al, 2004]. Furthermore, magnetic sepa¬ 
rators have been found to occur on a wide range of scales for 
they are found in global magnetic field extrapolations from 
photospheric synoptic magnetograms extending over several 
solar radii [e.g., Platten et al, 2014; Edwards et al, 2015] 
and also in local magnetic extrapolations of the quiet-Sun 
from high-resolution magnetograms [e.g.. Close et al, 2005]. 
Additionally, magnetic reconnection plays an important role 
in planetary magnetospheres enabling the interaction of the 
interplanetary and planetary magnetic fields and powering 
flux transfer events and substorms. Indeed, separators have 
been identified in numerous models of the Earth’s magne¬ 
tosphere [e.g., Hu et al, 2004; Laitinen et al, 2006, 2007; 
Dorelli et al, 2007; Dorelli and Bhattaeharjee, 2008, 2009; 
Hu et al, 2009; Pulkkinen et al, 2010; Ouellette et al, 2010; 
Peng et al, 2010; Cnossen et al, 2012; Komar et al, 2013]. 
Notably, Komar et al [2013] have shown that on the dayside 
magnetopause, regardless of the direction of the interplan¬ 
etary magnetic field (IMF), magnetic reconnection is most 
likely to occur at separators. 

There are only limited “observations” of magnetic sepa¬ 
rators in the Earth’s magnetosphere. This is not surprising 
since, as pointed out by Parnell et al [2010a], “for 3D mag¬ 
netic fields, global 3D topological structures and local 3D 
field structures do not [necessarily] coincide”. Since separa¬ 
tors are global structures this means it is difficult to deter¬ 
mine separators unless the magnetic field is known “practi¬ 
cally everywhere, which is obviously not the case from solar 
and magnetospheric observations”. Nonetheless published 
results claiming to have detected reconnection at separators 
using Cluster data include the work of Phan et al [2006] 
who reported an observed separator (called, in their paper, 
an X-line) on the Earth’s dayside, whilst Xiao et al [2007] 
report having observed separator reconnection in the night- 
side magnetosphere. Deng et al [2009] and Guo et al [2013] 
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also report examples of separator reconnection observed us¬ 
ing Cluster data. 

We note that in some of the works referenced here, the 
term ‘X-line’ is used in place of ‘magnetic separator’ since it 
is thought that the projected magnetic field in a plane per¬ 
pendicular to a 3D magnetic separator is X-type and looks 
like a 2D null point. This is not true since in 3D the global 
topology and local magnetic field are not necessarily coin¬ 
cident as they are in 2D. So although the projection of the 
global topology in a plane perpendicular to a separator is 
X-type, the projection of the local magnetic field in such a 
plane may be X-type or 0-type [Parnell et a/., 2010a]. Eur- 
ther to this, the term ‘X-line’ is often used to describe a line 
of nulls, which is an unstable feature [Hesse and Sehindler, 
1988] and entirely different from a magnetic separator. A 
3D magnetic separator may be regarded as a 2.5D X-line (a 
2D null point with a guide field). 

A separator is a special field line that connects a positive 
null to a negative null (separators may also connect two bald 
patches or a null and a bald patch [see, for example, Titov 
et al.^ 1993; Haynes, A. L., 2008], but we do not concern our¬ 
selves with such separators here). A positive/negative 3D 
null is associated with a pair of spine lines (field lines which 
are directed into/away from the null point) and a separa- 
trix surface (a surface of field lines pointing away from/into 
the null point). A separator is formed when; (i) the separa- 
trix surfaces of two oppositely signed nulls intersect, (ii) the 
spine of one null point intersects with the separatrix surface 
of another null point of the same sign and (iii) the spines 
of two oppositely signed nulls intersect, (i) is the only one 
which is both general and generic since the others are unsta¬ 
ble to perturbations (i.e., small perturbations would cause 
them to disappear resulting in a change to the structure of 
the topology of the magnetic field) and, hence, is the only 
type considered in this paper. A separator lies on the bound¬ 
ary of four topologically distinct flux domains (i.e., move¬ 
ment between the domains would result in a discontinuous 
jump in field line mapping) and so field lines in the vicinity 
of separators are sensitive to flows across these boundaries. 
Hence, currents build easily along separators [Lau and Finn, 
1990; Haynes et al, 2007; Parnell et al, 2010a, b; Steven¬ 
son et al, 2015] and so they are prime locations where 3D 
reconnection occurs. 

Separator reconnection has been studied both analyt¬ 
ically and numerically [Sonnerup, 1979; Lau and Finn, 
1990; Longeope and Cowley, 1996; Galsgaard and Nordlund, 
1997b; Galsgaard et al, 2000; Longeope, 2001; Pontin and 
Craig, 2006; Haynes et al, 2007; Dorelli and Bhattaeharjee, 
2008; Parnell et al, 2008, 2010a, b; Komar et al, 2013]. 
Reconnection at separators is known to be different from 
3D null point reconnection. In particular, Parnell et al 
[2010a] studied the dynamic nature of separator reconnec¬ 
tion in a model where two nulls, positioned on the base 
of the box, were moved together via boundary driving such 
that their separatrix surfaces intersected to form separators. 
They found that the parallel electric field along a separator 
varies spatially and temporally and may be multiply peaked, 
implying that there can be one or more local “hot spots” 
of reconnection along a separator. Additionally, these re¬ 
connection hot spots, which coincided with counter-rotating 
flows, occur away from the ends of the separator making 
separator reconnection distinct from 3D null point recon¬ 
nection. They discussed that the projected magnetic field 
in planes perpendicular to a separator may either be X-type 
or 0-type, with, in their model, the X-type projected field 
regions corresponding to 0-type flow and weaker reconnec¬ 
tion and the 0-type projected field regions corresponding 
to X-type flow and stronger reconnection. This finding has 
been questioned by some who believed the twisted nature 
of the field was a result of the specific driving in the exper¬ 
iment. 

In this work, we detail the properties of 3D separator re¬ 
connection in a non-driven experiment, to avoid any such 


problems, and also where the null points and separator are 
far from the boundary. The model starts from a MHS equi¬ 
librium, containing a separator current layer and excess en¬ 
ergy above that of a potential field, formed through the 
non-resistive relaxation of an initially non-potential, non¬ 
force-free field, discussed in detail in Stevenson et al [2015]. 
Reconnection is triggered at the separator current layer us¬ 
ing an anomalous diffusivity to mimic the onset of micro¬ 
instabilities. In nature it is believed that the slow driving 
of complex magnetic fields leads to equilibria forming which 
have current layers located, for instance, where the field line 
mapping is discontinuous. Reconnection can occur at such 
current layers via micro-instabilities once the length scales 
are sufficiently short such that the magnetic Reynolds num¬ 
ber, Rm < 1. 

This approach is different to many works, which start 
from potential (minimum energy) fields that are driven on 
the boundaries (at slow or fast rates) to induce reconnection. 
In these situations, the resulting reconnection rate is typi¬ 
cally found to depend on the rate of driving [e.g., Galsgaard 
and Parnell, 2005], and the nature of the initial magnetic 
configuration. 

The approach we use has been used before to investi¬ 
gate reconnection at 2D magnetic null points, e.g., Fuentes- 
Ferndndez et al [2012a, b] where high and low plasma-beta 
reconnection regimes were studied, respectively. The high 
and low-beta MHS configurations, from which the reconnec¬ 
tion experiments start, involve enhanced current not just at 
the null point, where it forms a null current layer, but also 
along the separ at rices of the null. Everywhere else the cur¬ 
rent is very low. In contrast to Longeope and Priest [2007] 
and Longeope and Tarr [2012] (zero-beta models), the nu¬ 
merical models of Puentes-Fernandez et al [2012a] (high- 
beta), and Puentes-Fernandez et al [2012b] (low-beta) find 
that the reconnection process converted most of the mag¬ 
netic energy (stored in the 2D null current layer configu¬ 
ration) directly into internal energy, via Ohmic dissipation, 
with only a little being converted initially into kinetic energy 
and then damped due to viscosity. Additionally, Fuentes- 
Ferndndez et al [2012a] found that the value of the mag¬ 
netic diffusivity affects not only the reconnection rate, but 
also the amount of magnetic energy converted into kinetic 
and internal energy. 

All experiments studied here start with equilibria in¬ 
volving a separator current layer embedded in a high-beta 
plasma (detailed in Stevenson et al [2015]). So called clus¬ 
ter separators (separators that link nulls clustered within 
a single weak field region [Parnell et al, 2010b]) are likely 
to be in a high-beta plasma. These separators are typi¬ 
cally short (approximately 1-2 Mm in length) [Parnell et al, 
2010b]. The vast majority of null points in the solar atmo¬ 
sphere occur low down [e.g., Regnier et al, 2008; Longeope 
and Parnell, 2009; Edwards and Parnell, 2015] and thus will 
reside in the high-beta chromospheric region of the Sun; in¬ 
tercluster separators (linking nulls from different null clus¬ 
ters) that lie within the chromosphere are likely to have 
lengths below 10 Mm (considering the typical size of small- 
scale photospheric magnetic features). The plasma beta in 
the magnetosphere lies between 1 and 10 [Trenehi et al, 
2008], thus separator reconnection here is also likely to oc¬ 
cur in a high-beta plasma. The lengths of these separators 
have been found to be greater than 200 Mm [Komar et al, 
2013; Komar, C. M., 2015]. As a result of this wide range 
of separator lengths, the results that we present here are 
non-dimensional. We explain the appropriate scalings that 
should be applied to produce dimensional results applicable 
to the many varied space physics situations where reconnec¬ 
tion occurs. 
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Here, we study the nature of non-driven 3D separa¬ 
tor reconnection at a generic single-separator current layer. 
Specifically, we choose an idealistic setup (a straight sin¬ 
gle separator) to allow us (i) to relatively easily analyse the 
nature of the reconnection found at all points along the sep¬ 
arator and (ii) the nature of the waves and flows [Stevenson 
and Parnell^ 2015] resulting from the reconnection. In par¬ 
ticular, we do not attempt to simulate one specific space 
physics event, but rather produce a model whose general 
results may be applicable in many events involving separa¬ 
tors. We look to answer: what is the partitioning of the 
magnetic energy released by the reconnection?; how is the 
reconnection rate affected by the value of the diffusivity, the 
size of the diffusion region and the background viscosity? 
We compare our results to those of Parnell et al. [2010a], 
e.g., is the reconnecting magnetic field elliptic (helical) or 
hyperbolic?; to what extent do the null points play a role in 
separator reconnection?; does reconnection occur along the 
entire length of the separator? 

In Sect. 2 we detail the properties of the MHS equilib¬ 
rium and the current layer it contains. The numerical model 
used to carry out the experiments is discussed in Sect. 3, 
followed by an analysis of the energetics of the main ex¬ 
periment (Sect. 4). We then examine the nature of the 
reconnection (Sect. 5) by analysing where the reconnection 
is strongest and the plasma properties there. Einally, we 
detail the effects of varying the diffusivity, the size of the 
diffusion region and the background viscosity in the system 
(Sect. 6) before summarising our results (Sect. 7). 

2. MHS Equilibrium Current Layer 

The MHS equilibrium used as the initial condition for 
our numerical experiments was formed through non-resistive 
MHD relaxation [Stevenson et aL, 2015]. The initial equi¬ 
librium is identical to that found by Stevenson et al. [2015] 
with the exception of an increased domain length in the 
^-direction, but the grid resolution is the same. The dimen¬ 
sions of the numerical domain, in Cartesian coordinates, are 
— 1.0 < x,y < 1.0 and —1.75 < < 2.75 and in grid cells is 

(512,512,768). 

The MHS equilibrium contains two 3D null points located 
at (0,0,-0.10) and (0,0,1.08) with a separator along the z-axis 
linking them created from the intersection of their separa- 
trix surfaces. (The null points in our model were tracked 
using the trilinear null finding method of Haynes and Par¬ 
nell [2007] and the magnetic skeleton was found using the 
method described in Haynes and Parnell [2010]). This equi¬ 
librium contains a twisted current layer which lies along the 
separator (Eig. la, a movie is available online showing a 
360° view of the magnetic skeleton with the current layer). 
The dominant component of current in the current layer is 
parallel to the separator. Eig. lb shows a horizontal cut 
perpendicular to the separator at z = 0.4 (the location of 
the peak current along the separator) in which the strong 
current at the separator (pink contour) and along the sepa- 
ratrix surfaces (cyan curved contours) is visible. In this cut, 
an insert of the current layer around the separator is in¬ 
cluded to indicate its width, le, and depth, d. We have also 


plotted a yellow contour at jcrit = 10: the current above 
which the diffusivity is non-zero (see Sect. 3). 
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Figure 1. (a) Skeleton of the MHS equilibrium magnetic 
field with purple isosurface of j|| = 10.0. Also shown are 
the positive/negative nulls (blue/red spheres) with asso¬ 
ciated spines (blue/red lines) and separatrix-surface field 
lines (pale-blue/pink lines) and the separator (green line, 
hidden by the current layer). The solid pale-blue/pink 
lines indicate where the separatrix surfaces intersect the 
boundaries. A movie is available online showing a 360° 
view, (b) Perpendicular cut across the MHS equilibrium 
separator at z = 0.4 showing contours of |j| and white 
solid/dashed lines going through the depth/across the 
width of the current layer, respectively. A yellow contour 
is drawn at jcrit = 10. The inserted image highlights the 
depth (d) and width (w) of the current layer in this plane. 

In the MHS equilibrium, the separatrix surfaces that cre¬ 
ate the separator are twisted about it forming cusp regions 
in planes perpendicular to the separator. Within these 
cusps the plasma pressure is enhanced and outwith them 
the plasma pressure is diminished. The equilibrium is not 
perfect: an infinite time would be required to achieve such a 
thing [Stevenson et a/., 2015]. Instead, very small localised 
residual pressure and Lorentz forces remain on the edge of 
the current layer about the separator and the separatrix 
surfaces near the separator. Outside the current layer the 
current is small and the field is in force balance. 
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3. Numerical Model 

The MHS equilibrium current layer has length Igep — 
1.18, depth d = 0.06 and width w = 0.24 (the latter two 
dimensions are calculated at z = 0.4). To study the recon¬ 
nection that can occur at this separator current layer due to 
micro-instabilities (modelled by introducing an anomalous 
diffusivity), we employ the 3D resistive MHD code LareSd 
[Arher et aL, 2001 ]. 

LareSd solves the MHD equations in a frame that moves 
with the fluid and then maps the resulting Lagrangian grid 
geometrically back onto the original Eulerian grid. R uses 
a staggered grid in which the pressure (p), internal energy 
per unit mass (e) and the density (p) are defined at the cell 
centres. The magnetic field components (Bx^By^Bz) are 
defined on the cell faces with the Evans and Hawley [1988] 
constrained transport method for magnetic flux employed 
to help maintain the solenoidal constraint, V • B = 0. The 
velocity components {vx^Vy^Vz) are defined on the cell ver¬ 
tices to prevent the checkerboard instability. The following 
normalised quantities (identified by the hats) represent di¬ 
mensionless variables used by the code. 


where ( 0 , 0 , ^z) and ( 0 , 0 , ^n) are the equilibrium positions of 
the lower and upper nulls respectively and Cs{z) and ca{z) 
are the sound speed and the Alfven speed at ( 0 , 0 ,z), re¬ 
spectively. 

We employ line-tied boundary conditions to prevent en¬ 
ergy leaving or entering the domain. Hence, the derivatives 
of the internal energy per unit mass and the density and all 
the components of the magnetic field normal to the bound¬ 
aries are set to zero. The velocity is also set to zero on the 
boundaries (v = 0 ). 

To gain reconnection, resistive terms are included in the 
governing equations. We choose to use a non-uniform diffu¬ 
sivity which is zero unless the current is greater than a set 
amount, jcrit, 

0 |j| ^ jcrit-) 

Vd |j|> jcrit. 


X = LnX, B = BnB and p = pnp, (1) 

where x = (x,p,z) is the position vector with Bn and 
Pn representing the normalised length, magnetic field and 
density, respectively. Using these normalising factors, the 
normalising constants for the velocity, pressure, current, in¬ 
ternal energy per unit mass and plasma beta may be written, 
respectively, as 


Vn 


B„ _ Bl 

, -7 Pn — 5 

y/P^Pn Mo 


jn — 



^n=Vn = and ^ ( 2 ) 

PoPn B^ 

Note, in the code’s dimensionless units the magnetic perme¬ 
ability Mo=l. This leads to the following resistive normalised 
MHD equations which are used in LareSd (note, the hats 
have been dropped for simplicity) 


Dp 

Dt 


-pV • V, 


(3) 


^ = i(V X B) X B - Np+ 
Dt p p p 


DB 


(B • V)v - B(V • v) - V X (pV X B), 


m 


■v+-H^ + 

P P P<^ 


(4) 

(5) 

( 6 ) 


where t is time, Yn — i^(V^v + |V(V • v)) is the vis¬ 
cous force (where v is the background viscosity), p = 
l/(Mocr) is the magnetic diffusivity where a is the electric 
conductivity, — ^(V • v)^) is the viscous 

heating term (where e^j, the rate of strain tensor, equals 
(dvi/dxj) + (dvj/dxi)) and a is the electric conductivity. 
The term j^/cr represents Ohmic dissipation and our clo¬ 
sure equation is e = pjp(pi — 1 ) where 7 , the ratio of specific 
heats, equals 5/3. The value of the background viscosity 
is V — 0.01. The background viscosity (which is equal to 
pvk where Vk is the kinematic viscosity) plays a role in the 
viscous force, iU, and viscous heating terms, discussed 
previously. 

All times in this paper are normalised with respect to t/, 
the time it would take a fast magnetoacoustic wave to travel 
from the lower null to the upper null along the path of the 
separator (z-axis) in the initial MHS equilibrium: 




f 

J Zi 




Csiz)^+c\{z) 


.dz = O.g 


(7) 


In our main experiment jcrit = 10 such that diffusion only 
occurs at the separator current layer (and not on the current 
on the separatrix surfaces), p^ = 0 . 001 , corresponding to an 
average Rm of 10^ along the separator, and u = 0.01. An 
analysis of how varying p^, jcrit and ly effect the reconnec¬ 
tion rate and the energy conversion is presented in Sect. 6 . 
The behaviour in all our experiments is very similar, so we 
mainly concentrate on one (the main) experiment in order to 
describe the basic behaviour that is seen in all experiments. 

The initial setup for all experiments discussed in this pa¬ 
per is a MHS equilibrium. When we apply a non-zero diffu¬ 
sivity, reconnection occurs immediately (as will be discussed 
in detail in Sect. 4). If p^ = 0, then the current layer in 
the equilibrium does not decay [see Stevenson et al.^ 2015]; 
indicating that the numerical diffusion is smaller than our 
range of values for pd- 

In the following sections, we describe the basic energetic 
behaviour and partitioning found in the main experiment 
(Sect. 4). Then the nature of the reconnection and magnetic 
field evolution are studied in Sect. 5, while Sect. 6 studies 
the effects of varying the diffusivity, the size of the diffusion 
region and the background viscosity on the energetics and 
reconnection rate. Einally, in Sect. 7, the conclusions are 
presented. 

4. Energetics 

In order to determine the response of the MHS equilib¬ 
rium to the introduction of an anomalous diffusivity we 
consider the system’s energetics. The change in energies 
against time for the main resistive 3D MHD experiment with 
Pd = 0.001, jcrit = 10 and v = 0.01 are plotted in Eig. 2a, 
where the magnetic, internal and kinetic energies are nor¬ 
malised to the maximum change in magnetic energy and the 
normalised kinetic energy is multiplied by fifty for represen¬ 
tational purposes. Note, time is plotted on a logarithmic 
scale. 
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Figure 2. Plots of (a) the change in energies, (b) the in¬ 
stantaneous heating terms and (c) the reconnection rate 
(blue line) and total flux reconnected (pink line). Note, 
the kinetic energy and viscous heating curves have been 
multiplied by fifty and both curves in (c) have been mul¬ 
tiplied by one hundred. The black dashed vertical lines 
highlight where the first phase ends (along with the sym¬ 
bols I and II) and the dashed horizontal line in (b) indi¬ 
cates where zero lies. 

Since there is no Poynting Elux on the boundaries (i.e., 
no gain or loss of energy through the boundaries) the change 
in total energy (AEt) should be zero throughout the exper¬ 
iment. We hnd that it is zero to within a relative error of 
2.1 X 10“^%, indicating that the code properly conserves 
energy. Eig. 2a shows that in the hrst 0.09t/ of the ex¬ 
periment, magnetic energy (AEb) is converted directly into 
internal energy (AEi) with only a little going into kinetic 
energy (AEk). During this time three quarters of the total 
loss in magnetic energy occurs. As seen from Eig. 2b, this 
magnetic energy is converted into internal energy via Ohmic 
dissipation with only 0.23% going into kinetic energy. 

Eig. 2c shows the reconnection rate and the total flux 
reconnected during the experiment (calculated from E\\dl 
along the separator as in Parnell et al. [2010a], an approach 
that is validated by Eig. 5). Clearly, the rate of reconnec¬ 
tion before t = 0.09t/ is dramatically different to that after 
and so we define this time as the change between two differ¬ 
ent phases: phase I, the fast-reconnection phase (88% of the 
hnal total reconnected flux is reconnected here) and phase 
II, the slow impulsive-bursty reconnection phase. Negligible 
reconnection occurs in the hrst few moments of the second 


phase (c.f. the hat gradient just after t = 0.09t/), but then 
small amounts of reconnection occur in bursts throughout 
the second phase. The vertical black dashed line on all these 
plots denotes where the change of phase occurs. 

The existence of two distinct reconnection phases is not 
unexpected since the 2D spontaneous reconnection exper¬ 
iments of Puentes-Fernandez et al. [2012a] and Fuentes- 
Ferndndez et al. [2012b] found the same behaviour. Fuentes- 
Ferndndez et al. [2012a] (high-beta) found that the fast- 
reconnection phase was followed by a steady, as opposed 
to the bursty, reconnection phase we hnd. In their low- 
beta model a bursty reconnection phase was found, suggest¬ 
ing that one of the key diherences between low and high 
beta spontaneous reconnection would be a much stronger 
implusive-bursty reconnection phase after the initial fast- 
reconnection phase. 

Diherences between the energetics of the two phases are 
evident in all graphs in Eig. 2. During phase I, in addition to 
the reconnection, which begins at a high rate, but decreases 
rapidly so it is small by the end of the hrst phase, there is a 
minimal contribution from the viscous heating term (whose 
peak during this phase is just (1/500)^^ of the peak Ohmic 
heating rate). A small amount of adiabatic cooling is also 
observed (Eig. 2b) indicating that locally a rapid expansion 
has occurred within the system. 

As the hrst phase (in which there is fast Ohmic dissipa¬ 
tion) comes to an end at t = 0.09t/, the gradients of the 
magnetic and internal energies change and the Ohmic heat¬ 
ing rate decreases: most of the current above jcrit in the 
separator current layer has been dissipated. In phase II, the 
kinetic energy, which until now has been slowly building, in¬ 
creases sharply associated with the presence of hows in the 
system following the rapid adiabatic cooling (see Stevenson 
and Parnell [2015] for further details of the waves and hows 
created by this reconnection). This increase in kinetic en¬ 
ergy is rehected by an increase in viscous heating to 8 x 10“^, 
however, the amount of Ohmic heating is always greater 
than the viscous heating throughout both phases due to the 
high value of the plasma beta (the mean value of the plasma 
beta is /3 = 4.8 at t = Ot/) which makes it hard for large 
waves to be produced. In Eig. 2b), the viscous heating has 
been multiplied by 50, which is why, in this plot, the vis¬ 
cous heating appears to be greater than the Ohmic heating 
in phase IT The kinetic energy, and viscous heating terms, 
level out towards the end of the experiment indicating that 
a phase of slow reconnection has started. 

The reconnection in this experiment, therefore, occurs 
in two phases: the hrst phase is a highly dynamic fast- 
reconnection phase with low velocities which is dominated 
by Ohmic heating (Ot/ < t < 0.09t/) and the second phase is 
a slow impulsive-bursty reconnection phase which has a com¬ 
paratively low level of Ohmic heating, but with higher ve¬ 
locities than in phase I providing greater amounts of viscous 
heating (0.09t/ < t < 0.76t/ where t = 0.76t/ is where we 
stop the experiment since at this time the waves, launched 
by the reconnection, approach the boundaries of the domain, 
see Stevenson and Parnell [2015] for more details). 

In this paper, we focus on the properties of the recon¬ 
nection in both phases of the experiments. In a follow up 
paper we will detail the properties of the waves which travel 
out from the diffusion site, due to the sudden lack of force 
balance, and set up flows in the system. 

5. Nature of the Reconnection 

5.1. Magnetic Field Evolution 

Reconnection at a 3D magnetic separator is found to 
have several similar characteristics to 2D null-point recon¬ 
nection (such as the flux from one pair of oppositely-situated 
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flux domains is transferred into another pair of oppositely- 
situated flux domains), shown in Eig. 3. Eig. 3 displays 
the nulls, their spines, the separator and the same sample 
field lines (drawn initially in the two flux domains outwith 
the cusp regions, i.e., Di and D 2 , respectively) at four times 
throughout the fast-reconnection phase, which lasts up un¬ 
til t = 0.09t/. The earliest two times are plotted from both 
side and top views to show which domains the field lines lie 
in before and after the reconnection occurs. In these plots, 
field lines which have not reconnected are drawn as black or 
grey lines, if they lie in domains Di and D 2 , respectively. 
Once reconnected, these field lines are coloured orange and 
grey and lie in domains D 3 and D 4 , respectively. This re¬ 
connection is not undertaken by a pairwise matching of held 
lines, but by continual reconnection of held lines during their 
passage through the diffusion region. 



with spines (blue/red lines) and dark green separator. 
Black and grey held lines are drawn, at t = Ot/, within 
the regions which lie outwith the cusps (i.e., Di and D 2 , 
respectively). When the grey/black held lines reconnect 
their colours change to orange/green and they lie in re¬ 
gions D 3 /D 4 , respectively, (a) and (c) t = Ot/ (side and 
top views), (b) and (d) t = 0.019t/ (side and top views), 
(e) t = 0.028t/ and (f) t = 0.09t/. 


As shown in Eig. 2, and conhrmed in Eig. 3, much of the 
reconnection occurs within the hrst phase. At t = 0.09t/, 
there are, however, some held lines which have not yet re¬ 
connected (e.g., the remaining black and grey held lines in 
Eig. 3f). These are an example of the held lines that are 
reconnected during the second phase of the reconnection. 
Movies, available online, show the evolution of all of these 
held lines from side and top views throughout the reconnec¬ 
tion experiment. 

One important point to note is that the black and grey 
(pre-reconnection) held lines do not lie along the separator, 
but cross by it, entering almost parallel to the spine of the 
lower null and then leaving parallel to the upper null’s spine. 
This means that the held lines twist through 180° about 
the separator before they reconnect. After they reconnect 
(orange and green lines) the held lines still run almost par¬ 
allel to the lower null’s spine, but now leave running parallel 
to the other side of the upper null’s spine. To do this the 


held lines run almost straight up the separator and barely 
twist around it at all. 

The reconnection at the separator rapidly dissipates the 
current above |j| = jcrit in the MHS equilibrium current 
layer, which lies along the length of the separator. The cur¬ 
rent at t = 0.09t/, in a cut across the separator at z = 0.4 
(Eig. 4, c.f. with Eig. lb which is the same plot, but in the 
initial equilibrium state) decreases at the separator current 
layer, but does not decrease along the separatrix surfaces 
since the value of current here is less than jcrit • The current 
along the entire length of the separator diminishes quickly 
from its peak of 22.4 and mean of 19.5 during the first phase 
until it is just below the value of jcrit at about |j| = 9.8 all 
the way along the separator. This is an indication that sep¬ 
arator reconnection occurs along nearly the whole length of 
the separator and not just at one point as in null-point re¬ 
connection (in agreement with Parnell et al. [2010a]). After 
the first phase, the current does not really change along the 
separator. 


t =0.09t 
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Figure 4. Contours of |j|, in a cut at 2 : = 0.4 across 
the separator, at t = 0.09t/. The intersections of the 
lower (pale-blue) and upper (pink) null’s separatrix sur¬ 
faces with this plane are also plotted at t = Ot/ (dashed 
lines) and at t = 0.09t/ (solid lines) on top of white lines 
so they are distinct from the contours. The separatrix 
surfaces split the cut up into four regions marked as Di, 
D 2 , D 3 and D 4 . The insert shows the separatrix surfaces 
of the equilibrium held and the held at t = 0.09t/, in the 
region 0.0 <x,y< 0 . 2 , close up. 

The dissipation of the current at the separator causes the 
separatrix surfaces of the nulls to open up slightly in the 
cusp regions, as hux is transferred from domains Di and D 2 
outside the cusps to the domains D 3 and D 4 that are inside 
the cusps. This behaviour is visible in Eig. 4 where the in¬ 
tersections of the separatrix surfaces of both nulls with this 
plane are plotted at t = Ot/ (dashed pale-blue/pink lines) 
and at the same time that the contours of |j| are drawn, 
t = 0.09t/ (solid pale-blue/pink lines). The separatrix sur¬ 
faces of the nulls are shown to open up slightly within the 
cusp regions (for closeup see Eig. 4 insert), as a consequence 
of the reconnection. 

The reconnection causes the nulls to move slightly apart 
in the ^-direction, with the maximum length of the separator 
throughout the experiment just 1.005 times greater than the 
equilibrium separator length. The lower null moves down¬ 
wards from its initial position, along the z-axis, by O.lLo 
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in the first 0.03t/, whereas the upper null moves upwards 
O.O 8 L 0 along the z-axis, in the first O.Olt/. After these times 
both nulls move very slowly away from each other (the lower 
null moves a further distance of 1 x 10 “^Lo and the upper 
null moves a further 5 x 10“^Lo) during the rest of the ex¬ 
periment. 

5.2. Nature of 

In 3D, a non-zero integral of Ey (= r]j\\) along a field 
line indicates the existence of reconnection [Schindler et a/., 
1988; Hornig and Schindler, 1996]. Eurthermore, in a situ¬ 
ation where there is a single simple diffusion region, such as 
we have here, the maximum integral of identifies the main 
reconnection site, as well as the reconnection rate [Hesse and 
Birn, 1993]. Since all the field lines that thread the non¬ 
ideal region {r]d ^ 0 ) thread planes that cross perpendicular 
to the separator, we determine the integral of Ey along field 
lines that thread the plane z = 0.4 at t = 0.019t/ during 
phase 1. A contour plot of J^Eyd/ (Eig. 5) shows that the 
strongest reconnection occurs at the separator (in agreement 
with Parnell et al. [2010a]). The insert in this figure shows a 
close up of the contours around the separator in this cut and 
highlights the strength of the localised peak of the integral 
of Ey along field lines at the separator. Weaker reconnection 
occurs on neighbouring field lines creating non-zero E\\dl 
regions along the separatrix surfaces that form the separa¬ 
tor. Now that we have confirmed that the reconnection in 
our experiment is separator reconnection, we look at the 
time evolution of various parameters along the separator to 
determine more details about the nature of the reconnection 
itself. 


t = 0.019tf 
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Figure 5. Contours of the integral of Ey along field 
lines which thread this plane perpendicular to the sepa¬ 
rator at z = 0.4, at t = 0.019t/. The insert shows a close 
up view of the region highlighted by the box around the 
separator. 

Eig. 6 a shows how the parallel component of the elec¬ 
tric field (Ey) along the separator evolves in time as the 
reconnection proceeds. The separator is normalised such 
that it has length one for all time according to the equation 
z* — [z — zi)llsep where zi (the ^-coordinate of the lower 
null) and Isep (the length of the separator) take their respec¬ 
tive values for the frame being considered. The strongest 
reconnection occurs between about z* = 0.16 and z* = 0.7. 
Midway between these points is close to where the initial 
peak in current lies indicating that the peak reconnection is 
occurring where the current is strongest, as expected. Note 
that the strongest values of Ey lie midway along the sepa¬ 
rator away from the nulls indicating that the null points are 


not involved in separator reconnection [in agreement with 
Parnell et al, 2010a]. 

After about t = 0.03t/, the peak Ey (Eig. 6 a) is about 
half of what it was at the start of the experiment and de¬ 
creases over the next t = 0.06t/ to almost nothing. During 
the second phase, patches of Ey >0 are visible in Eig. 6 a. 
These regions, which have |j| > jcrit, exist very briefly since 
the current in excess of jcrit is immediately dissipated due 
to the nature of the non-uniform diffusivity. Although these 
small reconnection events do not exist for long (compared to 
the reconnection in phase I) there are a significant number 
of these events such that the total flux reconnected slowly 
grows in phase II, as shown in Eig. 2c. 

The solid black curve in Eig. 6 a is the zeroth contour 
of the discriminant of the perpendicular component of the 
magnetic field (B_l), along the separator. This contour high¬ 
lights that the projected magnetic field lines are locally ellip¬ 
tic ( 0 -type), in planes perpendicular to the separator, in the 
regions where Ey is strongest, during phase 1. Elsewhere in 
phase I, and everywhere in phase II, the projected magnetic 
field lines are locally hyperbolic (X-type) in planes perpen¬ 
dicular to the separator. This behaviour is consistent with 
that discussed in Sec 5.1 where the pre-reconnected field 
wound around the separator 180°, whilst after reconnec¬ 
tion the field lines ran more parallel to the separator and 
did not cross it. This has been seen before in separator re¬ 
connection experiments by Parnell et aL [2010a] who found 
that cuts across the separator, where the reconnection was 
strongest, corresponded to 2D 0-type magnetic field config¬ 
urations. 

5.3. Behaviour of the Local Reconnection Flows 

Unlike in 2D reconnection, the flows local to separator re¬ 
connection are not simple stagnation flows. Schindler et al. 

1988]; Hesse and Schindler [1988]; Hornig and Schindler 

1996] have shown that, in 3D, the flow in the diffusion re¬ 
gion counter-rotates such that a continuum of field lines dis- 
continuously reconnect at the separator in a manner that is 
not pairwise. This is indeed the case for separator reconnec¬ 
tion, as evidenced by the vorticity parallel to the separator, 
= (V X v);s, (Fig. 6 b). Therefore, in 3D separator recon¬ 
nection the magnetic field lines undergo a rotational slippage 
about the separator, which is not surprising since they are 
effectively undoing their 180° twist about the separator. 

Indeed, Fig. 6 b shows that the streamlines curl around 
the separator in opposite directions above and below a 
point initially around z* = 0.48 (midway between where 
the current and pressure are initially highest on the sepa¬ 
rator). This point, about which the streamlines counter¬ 
rotate, moves downwards along the separator during phase 
I (when the value of Ey here is strong) before returning to 
around z* = 0.49 near the end of phase I where it remains 
during phase IT Thus, the point on the separator about 
which the flow counter-rotates is close to, but not coinci¬ 
dent with the position of the peak reconnection along the 
separator. 

To discover more about the nature of the flow we also look 
at the time evolution of the discriminant of the component 
of velocity perpendicular to the separator (v^) (coloured 
contours in Fig. 6 c). This plot indicates that, in planes 
perpendicular to the separator, the flow starts off 0 -type 
locally about the separator, although the velocities here are 
very small and so there is effectively no flow. The flow per¬ 
pendicular to the separator then becomes X-type away from 
the nulls during the first phase where strong reconnection 
is occurring (the black lines in Fig. 6 c are contours of Ey 
along the separator, as determined in Fig. 6 a). As the 
reconnection rate during phase I decreases the perpendicu¬ 
lar flow once again becomes 0-type briefly. For the second 
phase, the local perpendicular flow returns to X-type about 
the separator along almost its entire length except for very 
close to the upper null. 
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Figure 6. The time evolution of (a) E\\ (coloured contours) with the zeroth contour of the discriminant 
of over plotted as black lines and the nature of the 2D magnetic field lines in planes perpendicular 
to the separator annotated on the plot, (b) the component of the vorticity parallel to the separator 
{ujz = (V X v);s), (c) the discriminant of and contours of E\\ (shown in Eig. 6a) are over plotted as 
black lines and (d) the z-component of the velocity, along the separator (which is normalised to lie 
between z* = 0 and z* = 1). The vertical axis on each plot is logarithmic to highlight the behaviour 
during the first phase and the dashed line (and symbols I and II) indicate where the first phase ends and 
the second phase begins. 


Parnell et al. [2010a] studied the local velocity flow 
and local magnetic field structure in planes perpendic¬ 
ular to a 3D separator during separator reconnection 
and found a clear-cut relationship between velocity flow 
and magnetic field local to the separator: when the lo¬ 
cal flow is 0-type the local magnetic field is X-type, and 
vice versa. Here, comparing the contours of the discrim¬ 
inant of (Fig. 6c) with the contours of the discrimi¬ 
nant of Bx, black lines on Fig. 6a, we see that there is 
no clear correspondence between the nature of the local 
flow and field about the separator. At the start of the 
experiment, the magnetic field local to the separator is 
strongly 0-type (away from the nulls where, by defini¬ 
tion, it must be X-type), but the velocities are too small 
to really claim that the flow is actually 0-type. The O- 
type magnetic field persists on the separator and the ve¬ 


locity becomes X-type along almost the entire length of 
the separator (in agreement with Parnell et al. [2010a]). 
Briefly, when E\\ is decreasing rapidly to zero, the flow 
about the separator becomes 0-type: at this time the 
perpendicular magnetic field is already X-type (again in 
agreement with Parnell et al. [2010a]). After this an X- 
type velocity flow then returns and persists throughout 
the weak bursty reconnection phase (phase II) where 
the perpendicular magnetic field local to the separator 
is also X-type along its whole length (contrary to Par¬ 
nell et al. [2010a]). Thus, the relationship between the 
local flow and field within the diffusion region is not 
as simple as Parnell et al. [2010a] state, however, dur¬ 
ing rapid reconnection the configuration does tend to 
involve field lines with a twist that are untwisted by a 
counter-rotating stagnation-type flow. 
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We have considered the behaviour of the flows in to 
and out from the separator, but what are they like along 
the separator? (Fig. 6d) shows the ^-component of ve¬ 
locity along the separator, which, for the whole of phase 
I, is directed outwards towards the nulls from a point at 
z* = 0.63. This location is neither coincident with the 
location of the peak reconnection, nor with the point 
of counter-rotation, nor with the peak pressure along 
the separator (which is located at about z* = 0.57). 
The point of divergence of this flow abruptly changes 
as phase II starts and gradually moves up to the upper 
null such that along the separator the flow is almost 
entirely downwards. Towards the end of phase II, near 
the lower null an upflow along the separator is formed. 

So except for a brief period during which the recon¬ 
nection is rapidly decelerating towards the end of phase 
I, the local flow about the separator is essentially a form 
of counter-rotating stagnation flow in which the field in 
the diffusion region is swept into the separator whilst 
rotating and then outwards down towards the nulls and 
away from the separator, again whilst rotating. Due 
to the lack of co-location of the counter-rotation point 
and the outflow point of flow along the separator, in 
the region of peak reconnection the local flows are more 
complicated. Note, details of the nature of the flows fur¬ 
ther out from the separator are discussed in Stevenson 
and Parnell [2015]. 

Finally, the ^-component of the curl of B along the 
separator (not shown) decreases over time, revealing 
that the magnetic field relaxes as the reconnection pro¬ 
ceeds. So, the magnetic field behaviour is as expected 
and is consistent with the field relaxing, via reconnec¬ 
tion, from a twisted state. 

6. Parameter Analysis: jcrit and ly 

In order to investigate the effects of (i) the anoma¬ 
lous diffusivity rjd^ (ii) the size of the diffusion region (by 
changing the value of jcrit) and (hi) the background vis¬ 
cosity z/, we run additional experiments in which each 
of these parameters is varied in turn. The results from 
these runs are discussed in the subsections below. 

Experiments with three different values of diffusiv¬ 
ity {rjd = 0.0005, 0.001 - main experiment and 0.002) 
with jcrit = 10 and z/ = 0.01, then three different val¬ 
ues of jcrit {jcrit = 10 - main experiment, jcrit = 7.5 
and jcrit = 8-5) with rj = 0.001 and z/ = 0.01 and fi¬ 
nally three different values of z/ are studied (z/ = 0.005, 
0.01 -main experiment and 0 . 02 ) with 77 = 0.001 and 


jcrit = 10 are studied. None of the plasma parame¬ 
ters, nor the initial magnetic field have changed from 
the original experiment in these runs except for those 
stated here. 

6.1. Effects of Varying 77 ^, jcrit and zz on the 
Reconnection Rate 

First, we investigate how the reconnection rate and 
the total flux reconnected depend on r]d, jcrit and z/. As 
one might expect, the higher the value of r]d, the faster 
the reconnection, but the shorter the fast-reconnection 
phase (phase I) (Fig. 7a). The total amount of recon¬ 
nected flux in all three cases is similar, but appears to 
marginally increase as r]d increases, except in the experi¬ 
ment with rjd = 0.002. We believe this occurs since near 
the start of the experiment with the highest anomalous 
diffusivity {rjd = 0.002) additional nulls appear in the 
system. These nulls, which are formed in opposite-sign 
pairs, appear close to the locations of the original nulls 
and lead to the creation of extra intercluster separa¬ 
tors in addition to the original (intercluster) separator 
that linked the two original null points/null point clus¬ 
ters. To calculate the reconnection rate and total flux 
reconnected, we assume that there is just one reconnec¬ 
tion site; the original separator. However, when there 
are multiple separators/reconnection sites, this assump¬ 
tion breaks down, and so we may not be identifying all 
the reconnection that occurs in the experiment with 
r]d = 0.002. This is the only experiment discussed here 
which displays this behaviour. The creation of these 
multiple nulls and separators, and the reconnection as¬ 
sociated with them, will be studied in a follow up paper. 

From Fig. 7b, we see that the smaller the diffusion 
region (i.e., the higher the value of jcrit) the shorter the 
fast-reconnection phase (phase I). This, of course, is ex¬ 
pected as a smaller diffusion region contains less flux to 
reconnect. However, the peak reconnection rate is un¬ 
affected by the value of jcrit- Thus, the total amount 
of flux reconnected increases as jcrit is lowered. 

From Fig. 7c, we find that varying the value of z/ has 
little effect on the reconnection rate during the fast- 
reconnection phase. During phase H a smaller viscos¬ 
ity is associated with a marginally faster reconnection 
rate. This is because the resulting flows that drive the 
steady-state reconnection of phase H are stronger in a 
fluid that is less viscous. This leads to more flux being 
reconnected overall in the case with the lowest viscos¬ 
ity. 
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Figure 7. Eirst column: plots of the reconnection rate (solid lines) and the total flux reconnected 
(dashed lines), multiplied by 100, second column: the energy and third column: the heating terms all 
plotted against time for experiments with (top row) r]d = 0.0005, r]d = 0.001 and rjd = 0.002 with 
jcrit = 10 and u = 0.01, (middle row) jcrit = 7.5, jcrit = 8.5 and jcrit = 10 with r]d = 0.001 and ly = 0.01, 
(bottom row) ly = 0.005, ly = 0.01 and ly = 0.02 with r]d = 0.001 and jcrit = 10. In plots (a) to (c) dashed 
vertical lines (coloured to match the respective experiment) highlight the time at which phase I ends and 
phase II begins for each experiment. In plots (d) to (i) symbols (annotated on the plots) indicate where 
phase I ends and phase II begins. 


6.2. Effects of Varying Diffusivity rjd on the 
Energetics 

For the three experiments with varying the same 
basic behaviour of energies is found as that seen in the 
main experiment (Figs. 7d and 7g). The change in ener¬ 
gies are normalised to the maximum change in the mag¬ 
netic energy of the main experiment (where rjd = 0.001, 
jcrit = 10 and ly = 0.01) so that the energy release of all 
experiments can be compared. As before, most of the 
magnetic energy lost is converted directly into internal 
energy with only a little kinetic energy (multiplied by 
fifty for representational purposes) generated (Fig. 7d). 
Varying rjd leads to a change in the rate of loss of mag¬ 
netic energy: naturally, the experiment with the highest 
rjd experiences the most rapid loss, but the shortest fast- 


reconnection phase (phase I). The total loss of magnetic 
energy increases as rjd increases. 

Most of the additional energy comes directly from 
Ohmic heating. A small contribution comes from 
viscous heating (multiplied by fifty for representa¬ 
tional purposes), due to some wave damping occur¬ 
ring predominantly during phase II of each experiment 
(Fig. 7g). The rate of viscous heating is very simi¬ 
lar in all three experiments suggesting that the more 
rapid reconnection has not led to the creation of larger 
perturbations. In each experiment, a small amount of 
adiabatic cooling occurs in the later part of phase 1. 
This occurs due to the sudden expansion of the field 
about the separator current layer. 

From Fig. 7g, it appears that the peak Ohmic heat¬ 
ing rate is basically linearly proportional to rjd, as is the 
duration of phase 1. However, the duration of the main 
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Ohmic heating period (which, as shown in Sect. 4, is 
slightly longer than phase I) does not vary in the same 
way and is shorter for large r](i than expected from a 
linear fall off. The relative error in the total energy is 
small and of the same order for all three experiments. 

6.3. Effects of Varying jcrit on the Energetics 

Here, we consider the experiments with varying jcrit- 
From Fig. 7e, we see that having a lower value of jcrit 
means that the total loss in magnetic energy increases 
and it takes longer for the majority of the magnetic en¬ 
ergy to be converted into internal and kinetic energy. 
Furthermore, the kinetic energy is greater as jcrit de¬ 
creases. These differences are not surprising, since a 
lower jcrit creates a larger diffusion region in which more 
current can be dissipated and creates a greater loss in 
force balance of the system. We find that the free en¬ 
ergy released, (calculated as a percentage of the energy 
released if a uniform background diffusivity was used), 
increases as jcrit decreases, with 3% of the free energy 
released when jcrit = 10, 5% when jcrit = 8.5 and 8% 
when jcrit = 7.5. 

In Fig. 7h, it is no surprise to see that the initial 
amount of Ohmic heating increases as jcrit decreases. 
Additionally, a lower initial jcrit leads to greater vis¬ 
cous heating associated with the larger kinetic energies 
in the system. 

6.4. Effect of Varying ly on the Energetics 

The effects of varying the value of the background vis¬ 
cosity, z/, has little effect on the total loss in magnetic 
energy (Fig. 7f). This is reflected in the Ohmic dissipa¬ 
tion (Fig. 7i), which appears to be basically the same 
in all experiments. The proportion of magnetic energy 
converted to kinetic energy during the experiment does 
depend on viscosity with a larger z/, corresponding to 
greater viscous heating in phase I, leading to smaller 
kinetic energy since the waves are damped to a greater 
extent. This means that near the end of phase I, and 
throughout phase II, the kinetic energy is greatest for 
experiments with lower z/. In phase II, we find that the 
larger the z/ the weaker the Ohmic heating. The adia¬ 
batic cooling term appears to be unaffected by the size 
of the viscosity. 

7. Conclusions 

In this paper, we have studied the properties of spon¬ 
taneous (undriven) reconnection at a 3D separator cur¬ 
rent layer using a resistive MHD code. We start from a 
system containing free energy, which is in MHS equilib¬ 
rium everywhere save for very small forces at the current 
enhancements about the separator and separatrix sur¬ 
faces: a perfect equilibrium would take an infinite time 
to form [Stevenson et a/., 2015]. An anomalous diffusiv¬ 
ity, rjd^ is used, to mimic the onset of micro-instabilities, 
which only acts where the current is greater than a set 


amount, jcrit^ such that reconnection only occurs at the 
separator current layer and, for low enough jcrit ^ a little 
way along the separatrix surfaces. 

In our experiments, unlike driven separator reconnec¬ 
tion experiments, such as Haynes et al. [2007]; Parnell 
et al. [2010a], the reconnection starts immediately and 
its rate is dictated simply by the value of the anomalous 
diffusivity and is not influenced by an external driver. 
We find that the reconnection occurs in two distinct 
phases; a short fast-reconnection phase followed by a 
longer slow impulsive-bursty-reconnection phase. Such 
a partitioning into distinct phases is the same as that 
found in 2D spontaneous reconnection experiments at 
null points [e.g., Puentes-Fernandez et a/., 2012a, b], 
but here an impulsive bursty regime is found even at 
high beta. 

In the main experiment, where rjd = 0.001, jcrit = 10 
and z/ = 0.01, most (88%) of the reconnection occurs 
during the first short phase, that lasts just O.OOtj, in 
which the current is rapidly dissipated away from the 
current layer. During this phase magnetic energy is 
mainly (99.77%) converted directly into internal energy 
via Ohmic dissipation. Only a small amount (0.23%) 
of energy is converted from magnetic to kinetic energy 
during this phase and then transferred into internal en¬ 
ergy via viscous damping. Additionally, due to a rapid 
expansion of the plasma as a result of the sudden recon¬ 
nection, a small amount of adiabatic cooling is observed 
during phase 1. Even during the slow impulsive-bursty 
reconnection phase, the limited Ohmic heating domi¬ 
nates over the viscous heating. This is likely to be a 
consequence of the high-beta plasma within which the 
separator is embedded since, in such a plasma, as a 
current layer is formed a strong pressure gradient force 
rapidly develops balancing the Lorentz force making it 
harder for strong current layers to form than in the 
equivalent low-beta system. Consequently at a high- 
beta separator current layer the perturbation caused by 
the sudden onset of reconnection will be smaller than 
in the low-beta case. Furthermore, the combination of 
a smaller perturbation and denser plasma in the high- 
beta, as opposed to low-beta, system means that little 
plasma is accelerated and, hence, a low kinetic energy 
is found. Further studies are needed to establish un¬ 
der what conditions, if any, a low-beta regime produces 
greater viscous heating than Ohmic heating. 

Although the kinetic energy is relatively small, this 
in no way implies that waves and flows are not pro¬ 
duced. The rapid loss of equilibrium caused as a result 
of the reconnection in the separator current layer nat¬ 
urally launches waves that generate flows within the 
system. These waves and flows are discussed in detail 
in the second paper of this series Stevenson and Parnell 
[2015]. 

The rate of reconnection in the first phase is up to 
22 times faster than in the second phase, and the maxi¬ 
mum reconnection rate increases as rjd gets larger. How¬ 
ever, the duration of the first phase shortens as r]d in- 
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creases, such that the total amount of flux reconnected 
in this phase is independent of rjd- As jcrit increases, 
the length of the first phase also decreases, but here 
the total amount of flux reconnected in the first phase 
increases as jcrit decreases. Varying the value of the 
background viscosity has little effect on both the rate 
of reconnection and the length of the first phase. 

In this paper, the reconnection rates we quote are 
dimensionless. These rates may be scaled to dimen¬ 
sional values by multiplying by where Bn (T), 

Ln (m) and tn (s) are the normalising magnetic field, 
length and time. For example, numerical modelling of 
separator reconnection at the dayside magnetopause by 
Komar et al. [2013] and Komar, C. M. [2015] found sep¬ 
arators of length Ln = 223 Mm within magnetic fields 
with strengths Bn = 5 nT. Applying these values, with 
a time scale of = 1 s, we find the peak reconnection 
rate at the start of our main experiment is 5.6 x 10^ V. 
This value is over 100 times greater than the observed 
value found by Chisham et al. [2004] who looked at high 
latitude dayside reconnection with northward IMF. Our 
values better match with those found by Pinnock et al 
[2003] (which range from 1.1 x 10^ V to 9 x 10^ V) who 
found the rate of reconnection along the projection of a 
separator in the ionosphere. 

Similarly, our diffusivity {r]d = 0.001) can be scaled 
to dimensional values in the dayside magnetopause by 
multiplying by LnVA^/Rmi where vau is the Alfven 
speed scaling factor. Komar, C. M. [2015] found 
Alfven speeds of 380 kms“^, hence, our diffusivity 
r]d = LnVAn/Rm = 8.5 X 10^^ m^s“^ (where Rm = 10^ 
as discussed in Sect. 3). This value is comparable with 
that used in Komar et al [2013] {tjk = 6 x 10^^ m^s“^) 
who studied similar experiments to those detailed in 
Komar, C. M. [2015]. 

During phase I, the strength of the reconnection, 
which occurs asymmetrically along the entire length of 
the separator with the strongest reconnection occurring 
away from the null points, decreases and only weak, 
short-lived reconnection events occur in phase IT Dur¬ 
ing this second phase, the impulsive bursts of recon¬ 
nection occur randomly in small localised regions on 
the separator away from its ends. Each event is very 
weak and short-lived. These events do not start un¬ 
til t = O.lTt/, some 0.08t/ after the end of the first 
phase. During this time we suspect that the system is 
again trying to regain equilibrium and so builds up the 
current along the separator. The instant the current 
anywhere along the separator rises above jcrit a little 
burst of reconnection dissipates it. 

We have found that regions along the separator where 
the reconnection is strongest are associated with 0-type 
perpendicular magnetic field lines local to the separator 
and that regions where the reconnection is weaker, in 
phase I, and throughout all of phase II, are associated 
with perpendicular magnetic field that is locally X-type. 
Note that the signature of 0-type field actually comes 


from field lines that twist through just 180° as they run 
up the separator. On the other hand the X-type arises 
because the newly reconnected field lines slightly twist 
back on themselves as they run up the separator. There¬ 
fore, identifying a separator or separator reconnection 
from knowledge of the local magnetic field behaviour 
is not possible, as mentioned in the introduction: sep¬ 
arators are global features and therefore knowledge of 
the global field is required to locate them. However, 
since strong 3D reconnection is associated with a par¬ 
allel electric field and, hence, (under the conditions of a 
classical Ohm’s law) is associated with a parallel electric 
current, magnetic fields that locally have an 0-type na¬ 
ture in planes perpendicular to them are potential sites 
for reconnection. 

A counter-rotating flow also exists at a point along 
the separator close to, but not coincident with, the peak 
reconnection. This counter-rotation of the flow exists 
simultaneously with both 0-type and X-type flow lo¬ 
cal to the separator. Indeed, the strongest reconnection 
appears to be associated with a counter-rotating stag¬ 
nation flow about the separator with the reconnected 
field moving outwards towards the nulls and away from 
the separator. The system starts with very weak flow 
which appears initially to be elliptic perpendicular to 
the separator, but then turns to a stagnation-type flow, 
suggesting that the stagnation flow is generated as a 
result of the reconnection. 

In a companion paper [Stevenson and Parnell, 2015], 
we discuss the properties of the waves which are 
launched from the diffusion site as a consequence of the 
loss of force balance due to the reconnection. The flows 
which these waves set up are also discussed in detail. 
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